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Abstract. We consider the possibility that the Eddington-Born-Infeld (EBI) modified gravity 
provides an alternative explanation for the mass discrepancy in clusters of galaxies. For this 
purpose we derive the modified Einstein field equations, finding an additional ’’geometrical 
mass” term which provides an effective contribution to the gravitational binding energy. 
Using some approximations and assumptions for weak gravitational fields, and taking into 
account the collisionless relativistic Boltzmann equation, we derive a generalized version of 
the virial theorem in the framework of EBI gravity. We show that the ’’geometrical mass” 
term may account for the well known virial mass discrepancy in clusters of galaxies. We also 
derive the velocity dispersion relation for galaxies in the clusters, which could provide an 
efficient method for testing EBI gravity from astrophysical observations. 
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1 Introduction 

One of the major challenges of modern astrophysics is certainly to explain the mass grav¬ 
itationally bounded in structures like galaxies and clusters of galaxies. The incompatible 
behavior of the rotation curves of spiral galaxies with the theoretical prediction of New¬ 
ton’s gravity [1] - to which the Einstein’s theory reduces at this scale - as well as the mass 
discrepancy in clusters of galaxies [2-4], are usually explained by postulating the existence 
of a cold pressureless fluid, called dark matter, which interacts only gravitationally (see [5] 
for explanatory reviews on the dark matter properties). Despite of many attempts [6], the 
nature of dark matter is still unknown and the only convincing evidence for its existence is 
gravitational. In clusters of galaxies the strong bending of light indicates that there is a lot 
of matter heavily concentrated in a central region of the cluster. As is well known, the total 
mass of a cluster can be estimated at least in two ways. First, by considering the motions of 
the member galaxies of the cluster and using the virial theorem to estimate the virial mass 
My. Second, by adding the mass of each individual galaxy member of the cluster, the total 
baryonic mass is determined. It is found that My is much greater than M^,, with typical 
values of My/Mb ~ 20 — 30 [7-9]. This discrepancy is usually attributed to the existence 
of dark matter. Very important for gravitationally bounded systems, like the ones found in 
Astrophysics, the virial theorem is, however, dependent on the gravitational theory. Hence, 
an alternative way to attack the problem of discrepancy of the mass of clusters of galaxies is 
modifying the theory of gravity (for a thoroughly comprehensive review on modified gravity 
and its consequences see [10]). This issue has recently received a lot of attention and several 
versions of the virial theorem are found in the scientific literature by using modified gravity 
such as metric f{R) theory [11], Palatini f{R) theory [12], metric-Palatini hybrid gravity [13], 
brane-world models [14] and DGP-inspired L(R) gravity [15]. In this paper we study the 
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virial theorem in matter-gravity coupling modifications based on the so-called Eddington- 
Born-Infeld (EBI) gravity theory [16, 17] (for previous theoretical study and support about 
this theory see [18-20]). 

The EBI theory is based on Palatini variational formulation, the metric and connection 
being treated as independent variables. It is indeed a special class of bimetric theories 
of gravity^ and has already been used as an alternative to dark matter and dark energy, 
since the additional fields introduced can be interpreted as both components in different 
regimes [16, 22-25] (see Refs. [26] for functional extensions of EBI theory). EBI gravity has 
also been investigated in issues such as black hole geometries [27], structure and stability of 
compact stars [28, 29] and singularity-free cosmologies [17, 30], among many others. 

The aim of our work is to investigate if the geometric mass, arising from the extra 
component of the field equations, can explain the well-known virial theorem mass discrepancy 
in clusters of galaxies. Eor this purpose we use the collisionless Boltzmann equation in the 
EBI gravity field equations and derive a generalized virial theorem. We find that the total 
mass of a cluster of galaxies in the framework of EBI gravity is the sum of its baryonic mass 
and a geometric mass Mq which accounts for most of the mass of the cluster. 

The work plan is as follows: In Section 2, based on the Palatini variational approach, 
we present an outline of the generalized field equations in EBI theory, while in Section 3 we 
derive the basic equations for a static spherically symmetric gravitational field. The gener¬ 
alized virial theorem is deduced in Section 4, and in Section 5 we present some astrophysical 
applications. The results and conclusions of our work are presented in Section 6 . 


2 Eddington-Born-Infield gravity 


Let us consider the action coupled to an Eddington-Born-Infield theory as proposed by 
Bahados in [16]: 


S 


1 

167rG 




/ 


+ (rxLm{i',gf,y), 


( 2 . 1 ) 


where G is the Newton’s gravitational constant, R is the Ricci scalar, 7 is a dimensionless 
parameter and I is a scale. \A^^\ denotes the absolute value of the determinant for any tensor 
The metric g^i, couples to matter helds, and is the Ricci tensor constructed solely 
from the connections (7“^ and its derivatives: 

K,. = . ( 2 . 2 ) 

The action (2.1) is a functional of and of the independent connections which should 
not be confused with the Levi-Civita connections P^j, of the metric The matter La- 
grangian density depends only on the metric g^u and the matter helds ip. Varying the 
action ( 2 . 1 ) with respect to g^i, and we obtain the held equations 
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R/iu — ^ {g/ii' T iQ/ii'') 


(2.3) 

(2.4) 


'^See Refs. [21, 22] for discussions about the nature of bigravity theories such as stability, allowed parameter 
ranges and cosmological aspects. 
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where is a new metric satisfying the metricity condition Dp{yfqq^^'^) = 0, Dp is the 
covariant derivative built with the connections (7“^ and q is the determinant of q^y (see [16] 
for some details of calculation). The connections are then given by 

Cpy = {dyq^p + dpq^y - d^qpy). (2.5) 

Equation (2.3) is the modihed Einstein held equations, the second term on the right hand side 
being the contribution from the EBI action. In this framework the dynamic of the spacetime 
is described by equations (2.3) and (2.4). It is worth mentioning that in EBI gravity, one of 
the metric, g^y, couples to matter helds and has important physical meaning since the matter 
follow the geodesics given by connections built with g^y. Thus it is called physical metric 
or matter metric. The second metric, denoted by q^y, generates the symmetric connections 
Cpy which are used to build the tensor K^y. The connection C^y is the geometric connection 
because it determines the curvature of space-time, the matter fields are not coupled to the 
metric q^y (for studies on matter-gravity coupling and doubly coupling in a bimetric gravity 
see, e.g.. Refs. [31, 32]). 

3 Field equations for a system of identical and collisionless point particles 

Let us assume that the geometry of the cluster can be described by a time-oriented Lorentzian 
four-dimensional space-time manifold with spherical symmetry^. The metrics g^y and q^y of 
an isolated spherically symmetric cluster are given respectively by 

gpydx^dx'^ = —e'^^'^^dt^ + + r‘^d9‘^ + sin^ OdcfP'^ (3.1) 

qpydx^dx" = + r^d6^ + sin^ Odcj)^, (3.2) 

where ;z(r), A(r), r/(r) and a(r) are functions of the coordinate r only. The galaxies in 

the cluster are considered identical, collisionless point particles of mass m, and their space- 
time distribution is described by a distribution function fB{x^,u^), defined in the phase 
space, which obeys the general relativistic Boltzmann equation. The number of galaxies 
per unit volume n{x^,u^) of phase space, as well as the energy-momentum tensor of the 
matter Tpy, are determined by the distribution function /^(x^, u^) as n = f fsdu and Tpy = 
f fsixiUpUydu, respectively, where Up is the 4-velocity of each galaxy and du = durduQdu^/ut 
is the invariant volume element of the velocity space. Let (u?) {i = t, r, 0, </>) be the average 
value of the square of the components of the 4-velocity. This average is defined by {u^) = 
(1/n) f fBufdu. Then we have, for i t for instance, (u^) = (1/n) f fBU^du = Tu/mn. 
The mass density of galaxies is p = nm, so we have Tu = p{uf). Note that although we 
considered all the galaxies having the same mass m, the density p is not constant since it 
depends of the number density n{x^,u^) in the phase space. The components of the tensor 
Tpy are represented in terms of an effective density p^f / and an effective anisotropic pressure 
with components radial Pg// and tangential Pgjp which are defined as [34] 

Peff = p{u‘f), = = p{u‘l) = p{ul). (3.3) 

Taking into account this form of the energy-momentum tensor and the metrics given by (3.1) 
and (3.2), the gravitational held equations (2.3) describing a cluster of galaxies in EBI gravity 

^For a discussion on deviations from this assumption and the effects of a triaxial structure see [33]. 
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take the form 
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where ' = d/dr and " = d?/dr"^. By adding the gravitational field equations (3.4)-(3.7) we 
obtain the following equation 


/ v” v' v'\' 1 . , 2\ 
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where (tt^) = (w^) + (ti^) + (w^) + (tt^). The second group of equations (2.4) for the EBI 
gravity gives us 
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(3.10) 

(3.11) 

(3.12) 


Using these equations, the second term on the right hand side of equation (3.8) can be 
untangled supposing that p(r) and a(r) are slowly varying functions of the coordinate r (p' 
and a' small), such that we can neglect quadratic terms. Thus we rewrite (3.8) as 
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so that in the right hand side only the second metric (3.2) appears. 


(3.13) 


4 The virial theorem in Eddington-Born-Infield gravity 

Now we generalize the virial theorem to apply it to galaxy clusters which are described by the 
distribution function fs- This function, however, obeys the relativistic Boltzman differential 
equation which must be integrated over the velocity space and then, in conjunction with the 
gravitational field equation (3.13), provides the relativistic virial theorem for EBI gravity. 
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4.1 The relativistic Boltzmann equation 

The transport equation for the propagation of a collisionless system of particles in a curved 
Riemannian space-time is given by the relativistic Boltzmann equation [34, 35] 

= (4.1) 


where p" is the 4-momentum of the particle, and (i = 1 , 2 ,3) are the Christoffel symbols 
associated to the metric (3.1) which, in the EBI gravity, is the metric that couples to matter 
fields [25]. A simplification of the Boltzmann equation comes about if we introduce at any 
point X of the space-time an appropriate orthonormal frame of tetrads e“(x) (a = 0,1,2,3) 
satisfying the condition where is the Minkowski metric tensor. Any 

tangent vector at x can be expressed as = p°'ea, which defines the tetrad components 
p“. In the case of the spherically symmetric line element given by equation (3.1) we introduce 
the following frame of orthonormal vectors [34]: ej) = and 

= r sin 0(5^. Now, let be the 4-velocity of a typical galaxy, satisfying the condition 
= — 1, with tetrad components given by = u^e“. The relativistic Boltzmann equation 
(4.1) in tetrad components is given by 


a ,,9fB a b cdfB 

u -I- %^u u 
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dx>^ 


du^ 
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(4.2) 


where the distribution function fs = /b(x^,u“) and 7 ^^ = are the Ricci rotation 

coefficients [34]. Here we assume that the only coordinate dependence of the distribution 
function is upon the radial coordinate r, so that equation (4.2) becomes 
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cot 0 = 0 . 


(4.3) 


The spherical symmetry requires that the coefficient of cot 9 be zero, which implies that the 
distribution fs is a function of r,Ur and + only. Multiplying equation (4.3) by murdu, 
integrating over the velocity space and assuming that Jb vanishes sufficiently rapidly as the 
velocities tend to ± 00 , we obtain 

r) n Fiij 

^^[p('«r)] + f - P[{ul) + (ul) - 2{ul)] = 0. (4.4) 

Now we multiply equation (4.4) by 47 rr^, and integrate over the cluster radius to obtain 

1 r dv 

2Ar - - A7Tr^p{uf)—dr = 0, (4.5) 

where 

K = f 2 -kp [{uD + {uD + {ul)]r‘^dr (4.6) 

Jo 

is the kinetic energy of the cluster, which is a fundamental piece in the virial theorem. 
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4.2 Geometric quantities 

In what follows we introduce some approximations that apply to test particles in stable 
motion around central fields and define our geometric quantities based on the q^i, metric. 
Firstly, we assume that u{r) and A(r) are small and slowly varying functions of the radial 
coordinate, so that in the left hand side of equation (3.13) the quadratic derivative terms 
can be neglected. Secondly, we assume that the galaxies in the clusters have velocities much 
smaller than the velocity of the light, so that (u^) « {uq) ^ (u^) <C (u^) « 1. Thus, equations 

(3.13) and (4.5) become, respectively 


= 4:7^0 p + AirGpg, 


(4.7) 



(4.8) 


where we have defined a geometric density pq as 



(4.9) 


Comparing (4.7) and (4.9) we see that the second metric (3.2) contribution appears as an 


additional source term. The total mass of the system inside a radius r is given by 


M(r) = / Airpr'^dr. 

Jo 

Therefore, multiplying equation (4.7) by and integrating from 0 to r we obtain 


(4.10) 



(4.11) 


where Mg is the geometric mass of the cluster, inside a radius r, defined by 



(4.12) 


We associate to this geometric mass a geometric potential energy defined as 



(4.13) 


and a geometric radius Rg of the cluster 



(4.14) 


In terms of this quantities we have 



(4.15) 


a relation which will be useful in what follows. 
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4.3 The virial theorem 

Finally we multiply equation (4.11) by dM{r)/r and integrating from 0 to R, we obtain 


2ii: + O + Og = 0, 


(4.16) 


which is the generalized virial theorem in FBI gravity, fig is given by (4.13), 14 is the 
gravitational potential energy due to the baryonic mass, given by 


14 = - 



GM{r) 

r 


dM{r ), 


(4.17) 


and we used equation (4.8) for the kinetic energy. In order to translate the virial theorem 
(4.16) into an expression for M^/M, we introduce the virial radius R^, defined as [34] 


Rv = 




J^MndMir) 

such that 14 = —GM'^/R^. In addition, taking the virial mass Mv, defined as [34] 


(4.18) 


2K 


GM^ 

Rv 


(4.19) 


and substituting these definitions in (4.16), taking into account (4.15), 
theorem as _ 


Mv _ 

M Y ^ R^\M J 


we rewrite the virial 


(4.20) 


For most of the observed galactic clusters, the relation Mv/M > 3 is true. Therefore, the 
unity term can be neglected in (4.20) and the virial mass in FBI gravity can be approximated 
by 


Mv 



(4.21) 


According to equation (4.21), most of the mass in a cluster with mass Mtot is in the form 
of the geometric mass Mg, so that Mq ~ Mtot- In other words, the ratio of the total mass 
and of baryonic mass is determined by a purely geometric quantity. The gravitational effects 
associated to the presence of the metric could in principle be tested through gravitational 
leasing. 


5 Astrophysical applications 

5.1 Estimating the geometric mass for galactic clusters 

Clusters of galaxies are the largest and most massive self-gravitating bounded systems in 
the Universe. Although clusters are dynamically evolving, and deviations from hydrostatic 
and virial equilibrium must be expected, simulations from Gaussian random density fields 
predict surprisingly tight virial relation [36] (see also [37]). Here we intend to estimate the 
order of magnitude of the geometric mass, hence we choose the simplest relations in modeling 
clusters. 
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It is well known that inside clusters of galaxies there is a hot tenuous gas, observed by 
its emission in X-rays predominantly through thermal Bremsstrahlung [38]. The following 
equation for the intracluster gas density pg, known as /3-model, provides a reasonably good 
description of the observational data [39]: 

/ 2 \ -3/3/2 

Pgir) = /^o (^1 + y (5.1) 

where Vc is the core radius, po and [3 are cluster-dependent constants. The mass Mg{r) of 
the gas inside a radius r is then given by 

r r r'^dr 

Mgir) = A^pgr^dr = dvrpo ^ ^2/^2)3/3/2 ' (^'2) 


The total mass Mtotir) inside the cluster can be obtained as function of the gas density via 
Jeans’s equation for a spherical system [9] 


Pg dr 


{pg{ul)) + ^ {2{ul) - {uf) - (ul)) 


d^ 
dr ’ 


(5.3) 


where <I>(r) is the gravitational potential. Assuming that the gas in the cluster is isotropically 
distributed, we take = (ttg) = (u^), and the pressure Pg = pg{u^). Since the gravitational 
field inside the cluster is weak, we assume that the gravitational potential satisfies the Poisson 
equation V^$(r) = dvrGptot) where ptot includes energy density of other forms of matter 
different from gas, like luminous matter, massive neutrinos, etc. Besides, the observed X- 
ray emission from the intracluster gas is usually interpreted by assuming that the gas is 
in isothermal equilibrium [38, 39]. Therefore, one may assume that the pressure Pg of the 
gas satisfies the ideal gas equation of state Pg = {k^Tg/pmpjpg, where is Boltzmann’s 
constant, Tg is the gas temperature, p ~ 0.61 is the mean atomic weight of the particles 
in the gas, and nip is the proton mass. A first integration of the Poisson equation gives 
d^/dr = CMtotif)/r'^. Using (5.1) and (5.3) we obtain the total mass profile inside radius r 


as 


Mtot{r) 


SksPPg r 
pnipG iPr'l/r'^ 


(5.4) 


On the other hand, the total mass of the cluster, according to the modified EBI gravity, 
consists of the sum of the baryonic mass (mainly the intra-cluster gas), and the geometric 
mass, so that 


pr 

Mtot(r) = 47r / {pg + pq)r^dr. (5.5) 

Jo 

It follows that the total mass inside the radius r satisfies the following mass continuity 
equation 

= 47rr^pg(r) -|- 47rr^pg(r). (5.6) 

Therefore we obtain the expression for the geometric density inside the cluster by using 
equations (5.1) and (5.4) 


_ 3ksf3Tg (l-b3rg/r^) 1 _ Po 

AirCpnip (1 -|- r^/r^)^ (1 -|- jp-ppPP 


(5.7) 
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as well as the geometric mass: 


M,(r) 


3kB(3Tg r A r 

Gfiirip l + r^/r^ Jo (1 + r2/r2)3/3/2 


(5.8) 


In order to estimate the geometric mass Mq{r), we consider regions where r S> Tc- In this 
approach the equations (5.7) and (5.8) reduces to 


Pqir) 


Sk^/^Tg I 

AirGiimp 


and 


T\,T ( \ (^^B^Tg 


r — 


(rc\^P 

- ( 7 ] 

(5.9) 

47rporf \ oa-0) 

3(1-/3)i 

(5.10) 


respectively. If we neglect the contribution of the gas in comparing with the gravitational 
effect of the geometric mass due to EBI gravity, then we approximate the above two equations 
by 


and 


Pq{r) « 

/ Sk^/STg \ 1 

\47rGpmp J 

(5.11) 

Mq{r) : 

f3kBf3Tg\ 

V pGrUp J ^ ■ 

(5.12) 


In this case, an upper bound for the cutoff of Mq{r) may be estimated if we consider the 
coordinate radius for which the decaying density profile (5.11) becomes equal to the mean 
energy density of the Universe puniv Let us name this radius . So, assuming pq{Rq) = 
Puniv = SR/SttG = 4.6975 x lO^^^/igQ g/cm^ [38] where H = bOh^o km/Mpc/s, and noting 
that k^Tg « 5 keV for most clusters, we obtain 


r>{cr) 

± Lg 


SjSksTg 


Y iTrGpmpPnniv 
= 25.06\/B \/ f hrnMpc. 


(5.13) 


By using (5.13) we find that the total geometric mass corresponding to this radius is 

/ h T \ 

= 13.7 X 101''^3/2 hg-o'Mo. (5.14) 

Considering /3 Ri 2/3 and k^Tg = 5 keV [38] we obtain = 20.46 h^g^Mpc = 15Mpc and 
Mq^''^ = 7.5 X 10^^/i^q^Mq = 5.5 x IO^^Mq for = 50/67.8, a value which is consistent 
with observations [38]. 


5.2 Typical values for the virial mass and radius 

Astrophysical observations, together with cosmological simulations, are generally interpreted 
in terms of a virialized part of the cluster. The radii commonly used are either i ?200 or 
i? 50 o- These radii corresponds to a fixed density contrast 5 « 200 (or 5 ^ 500) as compared 
to the critical density of the universe Puniv The corresponding masses inside this radii are 
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defined as M 200 and M 500 and it is usually assumed that the virial mass of the cluster is 
My = M 200 (or My = M 500 ) and the virial radius Ry = -R 200 (or Ry = -R 500 ) [38]. In order 
to compare the predictions of EBI gravity with the observations we estimate the virialized 
mass My{r = Ry) using (5.9) 

5 X X Rl^, = - Po , (5.15) 

and taking into account that for r ^ rc (5.1) give us poi^c/R 20 o)^^ = Pg{R 20 o)- We obtain 


R 200 


3/3kBTg 


47rG/i?Tip/3univ(<^ “1“ <^g) 


per 


(5.16) 


where 6g = Pg{R 20 o)/Puniv gives the density contrast of the intracluster gas as compared to 
the critical density of the universe at the radius i? 200 - Using the value found above for R^^ 
and taking 6g = 20, we obtain 

R 200 = 1.38/i^o^Mpc, ^ M 200 = 5.0 X lO^^/i^o^M©. (5.17) 

This value for 72200 is comparable to the analogous value R^^ = 3.51/i^Q^Mpc found in [14] 
for brane world models. Taking h'^Q = 50/67.8, Eq. (5.17) give us R 200 — 1.0 Mpe and 
Af2oo — 3.6 X lO^'^M©. A similar calculation, using now 6g = PgiR^oo)/Puniv = 50, give us 

72500 = 0.87/i5o^Mpc, ^ M 500 = 3.2 x IO^^/i^q^M©. (5.18) 

Observations show that the mass in the clusters of galaxies range from IO^^/i^q^M© to 
lO^^/igg^M© [38, 40], therefore the value of the geometric mass obtained in the framework of 
EBI gravity is consistent with observations^. 


5.3 Radial velocity dispersion in galactic clusters 

From the observational viewpoint the virial mass My is determined from the study of the 
velocity dispersion of the stars and of the galaxies in the clusters. In terms of the velocity 
dispersion ai the virial mass can also be expressed as [41] 

My = ^Ry. (5.19) 

Ct 

As done earlier, let us assume that the velocity distribution in the cluster is isotropic, that 
is: (u^) = (ui) + (u^) + (ug) = 3(uf) = 3(T^, where is the radial velocity dispersion and is 
related to ai by = 3af. Under this assumption, the radial velocity dispersion relation for 

clusters of galaxies in EBI gravity can be derived from equation (4.4), rewritten as 

-^{paj) + ^pu' = 0. (5.20) 

A first integration of the equation (4.7) yields 

= GMq{r) + GM{r) + G, (5.21) 

^The reader should be aware that in some of this references the virial masses are displayed in terms of 
^70 (= 1 - 4^50 )■ 
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where C is an arbitrary constant of integration. Substituting (5.21) into (5.20) we obtain the 
solution 

al{r) = -- r [GMg{r) + GM{r) + C] (5.22) 

P Jo 

As an example, let us consider a simple case in which the density p of the matter inside the 
cluster has a power law distribution given by 


p{r) = pqt 


(5.23) 


where po and n are positive constants. The corresponding mass prohle is M (r) = A.'KpQr^~'^ j (3— 
n). Assuming for the geometric mass GMq{r) = q^r where qo = ^ks^Tg/pnip, we obtain the 
following expressions for the velocity dispersion: 



C Cl 

qo + - -27rGpo r- In r H-r, 

2r po 


(5.24) 


for n = 1; 


for n = 3; and 



qo 
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ttGpo 


^Inr + 



1 Cl Cl , 

~I “I- ~A -1- 

4 r Po 


2/ 1 _ 90 , 27rCpo 

n^(n-l)(3-n) 


^2-n^ 


Cl Cl , 

-h —C 

n + 1 r Po 


(5.25) 

(5.26) 


for n 7 ^ 1,3, where we have introduced a new constant of integration Ci. The equations 
(5.24)-(5.26) are the velocity dispersion profiles predicted by EBI gravity and it can be used 
to estimate the virial mass in clusters of galaxies. It should be noted that in equations 
(5.24)-(5.26) the velocity dispersion shows the same radial dependency as that obtained in 
the framework of metric f{R) gravity by [11] and for brane world models by [14]. Therefore, 
measurements of radial velocity dispersions alone can not distinguish these three gravity 
models. On the other hand, the corresponding equations obtained by [12] in the Palatini 
formalism for f{R) gravity, as well as those obtained by [13] and [15] for hybrid metric- 
Palatini gravity and warped DGP-inspired L{R) gravity respectively, have mismatch in at 
least one term. This can be used, together with observational data, to test the different 
predictions of these theories. In general, the observational data are htted with these functions 
by using a nonlinear fitting procedure (see e.g. [41]). 


6 Discussions and final remarks 

Dark matter is a fundamental ingredient of the modern Cosmology, without which it seems 
impossible to explain the formation of structures in the Universe. However, although many 
efforts have been made so far, there is no observational evidence of non-gravitational inter¬ 
actions for this type of matter; accelerator and decaying experiments give no support for the 
physics upon which the dark matter hypothesis is based. This, indeed, rises doubts about its 
physical existence and, in turn, opens the possibility that the Einstein’s theory of gravitation 
breaks down at some scale (see [42] for development of quantitative procedures for comparing 
tests of modern gravity theories on all scales). 

The recent proposal by Bahados et al. [16, 17, 25], based on the Eddington-inspired- 
Born-Infeld gravity [18-20], has attracted considerable interest lately as a theory capable of 
explaining some of the dark sector phenomena. In particular, this theory can challenge the 


- 11 - 





necessity of including huge amounts of dark matter to explain the mass discrepancy in clusters 
of galaxies. In the present paper we analyzed this problem in the context of EBI gravity 
by deriving the untangled modified gravitational field equations (3.13). The extra terms, 
appearing in the modified gravitational field equations, induce an additional gravitational 
interaction, which in principle can account for the missing mass in clusters of galaxies. Using 
reasonable approximations for weak central fields, and taking into account the collisionless 
Boltzmann equation, we have derived a generalized version of the virial theorem within the 
context of EBI gravity (see Eq. (4.16)). The new virial mass is mainly determined by the 
geometric mass associated with the geometrical terms (see Eqs. (4.20) and (4.21)), showing 
the existence of a relationship of proportionality between the virial mass and geometrical 
mass of the cluster. Using the simplest relations and assumptions in modeling clusters of 
galaxies, as well as the intra-cluster hot gas, we have estimated the order of magnitude of 
the geometric mass, showing that it is compatible with observations (see Eqs. (5.17) and 
(5.18)). In order to compare our results with the approximations made in section (5) we write 

Eq. (4.21) for Rq = R!'q'^\ giving the relation /My = /RyJ ■ In subsection 

(5.1) we estimated a critical radius R^^'^ taking into account the mean energy density of the 
Universe puniv today, obtaining Rq^^'^ = 15 Mpc. Thus, taking Ry = R 200 = 1 Mpc, we obtain 
/M 200 ~ 3.87, a reasonable value since is un upper bound for the cutoff of Mq. 

On the other hand, as a result of our approximations in subsection (5.2) we obtained that 
R^q''^ can also be given by R^'q'^^ = {6 + 6 gY/‘^R 2 m = (200 + 20)^/^i720o (see Eq. (5.16)). Using 
again Eq. (4.21) with this value of R^q'^'^ we obtain now / M 200 ~ 3.85, which is almost 

the same value calculated above. A similar calculation, now with M 500 (see Eqs. (5.18)), give 
us ~ 4:.84 in the two approximations (remember now R^g'"^ = (500 + 50)^/^ii5oo)- 

The similarity of this values shows the consistency of our result (4.21) as well as the adequacy 
of the approximations made in section (5), specially in subsection (5.2). 

Our study of the virial theorem, generalized in the context of EBI gravity, makes evident 
the existence of a geometric mass (gravitational interaction) which would not be excluded 
by astrophysical observations. This might be an efficient tool to test the viability of this 
class of modihed gravity. Einally, recalling that the virial mass My is obtained from the 
observational study of the velocity dispersions of the stars in the cluster, we have derived a 
relation for the radial velocity dispersion (5.22) in the context of EBI gravity, which can be 
used for estimating virial masses. 
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